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Abstract 



The nonnegativity of the density operator of a state is faithfuhy coded in its 
Wigner distribution, and this places constraints on the moments of the Wigner 
distribution. These constraints are presented in a canonically invariant form 
which is both concise and exphcit. Since the conventional uncertainty princi- 
ple is such a constraint on the first and second moments, our result constitutes 
a generalization of the same to all orders. Possible application in quantum 
state reconstruction using optical homodyne tomography is noted. 
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The uncertainty principle exhibits a fundamental manner in which the quantum descrip- 
tion of nature departs from the classical one. For the canonical pair of variables {q,p) the 
Heisenberg commutation relation [q,p] = ih leads, for any state to the unbeatable 
limitation 

where (g) = {ip \q\ ip), Aq = q — (q), and so on. Every Gaussian pure state saturates this 
inequality. An important attribute of the uncertainty principle (|l]) is that it is invariant 
under all real linear canonical transformations, just as the canonical commutation relation 
is. 

This inequality can be generalized in a naive manner to higher orders in q, p. For any 
pair of hermitian operators A, B and state {ip) we have the Schwartz inequality 




(2) 

It is saturated if an only if Ali/j) and B \ip) are linearly dependent as vectors. Clearly, (|l]) 
is a particular case of (||) corresponding to A = q — {q) , B = p — (p) . Clearly, the choice 
A = q"^ — {q"^), B = p^ — (p^) will lead to a higher order uncertainty principle involving (g"^), 
{p'^)] the Fock states In), being eigenstates of (f +p^, will be expected to saturate this higher 
order uncertainty principle. That they indeed do so can be explicitly verified. 

One may indeed produce any number of such naive generalized uncertainty principles by 
making various choices for A, i? in (|^). But every one of them will suffer from the deficiency 
of not being invariant under linear canonical transformations. Further, there seems to be no 
reasonable sense in which the set of all such generalizations based on (j^) can be considered 
to be complete. 

The purpose of this Letter is to present a generalization of the uncertainty principle 
which largely overcomes these difficulties. This is achieved by applying to the Wigner 
quasiprobability concepts and results from the classical problem of moments [^. The 
final result is a nested sequence of constraints on the moments of the Wigner distribution. 
These constraints are tailored to capture the positivity of the density operator of a quantum 
state. Equivalently, a given real phase space distribution has to necessarily meet these 
constraints in order to qualify to be a bonafide Wigner distribution. 

It should be appreciated that the higher moments of the Wigner distribution are no more 
objects of purely academic interest. An enormous progress in quantum state reconstruction 
using optical homodyde tomography has been achieved in the last few years: the Wigner 
distribution of a state can now be fully mapped out , as has been demonstrated by several 
groups 0-0. 

There exist rigorous and mathematically sophisticated approaches to the quantum me- 
chanical moment problem But our considerations here are explicit and take full advan- 
tage of the canonical invariance underlying the Heisenberg commutation relation. 

Details of a classical probability density p{x) are coded in its moments 7n = / dxx^p^x). 
An important result in the problem of moments is this 0: given a sequence of numbers it 
qualifies to be the moment sequence of a bonafide probability distribution if and only if the 
symmetric matrix defined below is nonnegative: 
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>0. (3) 



This can be broken into a sequence of positivity conditions on the determinants of the 
submatrices of F, which in turn can be viewed as a nested sequence of constraints on the 
moments 7„; and these constraints are tailored to capture the pointwise nonnegativity of 
p{x). Reconstruction of p{x) from its moment sequence is the other part of the classical 
problem of moments 0. 

In quantum mechanics, the state is described not by a true probability density in phase 
space, but by one of several possible quasiprobabilities The earliest, and probably the 
most prominent, quasiprobability is the one introduced by Wigner It is intimately related 
to the Weyl ordering rule of association between the algebra A of functions / {q,p) of the 
phase space variables and the algebra A of operator valued functions F (g, p) of the canonical 
operators. The rule is specified first through the one to one correspondence e^'^^'^^ < — > e^^+'^P 
for plane waves, and then extended linearly to the entire algebra using Fourier techniques. 

The Weyl rule could equally well be specified in the monomial basis instead of the plane 
wave basis through the association q"^p^ < — > Tm,n for m, n = 0, 1, 2, ■ ■ ■ where the Weyl 
ordered monomial Tm,n is the coefficient of (m!n!)~^ ffm^n Taylor expansion of e^^^'^^. 

This is an isomorphism between A and A only at the level of vector spaces but not at the 
level of algebras. In particular, the product of two T^.n's is not another monomial but a 
linear combination of monomials 0]: 

Tm,nTm' ,n' ^ ^ -^m+m'— r— s ,n+n'— r— si 

s+r 



-ir(f)' 



'^'•'^ {m-sy.{n-r)\ [ r )[s ) ' 
The intimate connection between Weyl ordering and Wigner distribution is this: 

tT{pfm,n) = J dqdpq^'p^Wiq^p) . (5) 

That is, the quantum mechanical expectation of the Weyl ordered monomial Tm,n is precisely 
the mn-th moment of the Wigner function. By linearity, similar relation holds for any pair 
/ {q,p), F {q,p) related by Weyl ordering. 

The monomials Tm^n are hermitian, and transform in a simple manner under the group 
Sp (2, 3?) of real linear canonical transformations. This group can be identified with SL (2, 3f?), 
the group of 2 x 2 real matrices with unit determinant. Sp (2, 3?) acts identically on the pairs 
{q,p) and {q,p), and this action induces linear transformation in the algebras A and A in 
the natural manner. 

The set of homogeneous polynomials of order 2j in q and p (being linear combinations 
of qi~^pi^'^ for s = —j, + transform linearly among themselves under this 

transformation, leading to the spin-j representation of Sp(2,3?) in A. The Tm,nS in A 
transform in the same manner as the g^p^'s in A, and thus the vector space A decouples 
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into a direct sum of invariant subspaces under Sp(2,3fJ): A = V^^'^ © V^^'^ © V^^'^ © ••• 
Clearly, V^^^ is of dimension 2j + 1, and is spanned by = Tj_.s,j+s with s running over 
the range s = —j, + It acts as the carrier space for the spin-j representation 

of Sp (2, 3?) in A. Thus, every spin-j representation of Sp (2, 3?) occurs in A once and only 
once. 

It is convenient to arrange the ^j,s's for fixed j into a 2j + 1 dimensional column vector 

- (j) 

^ and then, for any chosen J, arrange these columns into a grand column vector of 
dimension (J + 1)(2J + 1). 

Let the {2j + 1) x (2j + 1) matrix K^^'^ (S) denote the spin-j representation for S G 
Sp (2, 3?). Since the defining representation of Sp (2, 3?) is the spin-| representation, we have 
K^^'^ (S) = S. Let Kj (S) be the block diagonal matrix of order ( J + 1) x (2 J + 1) with 
diagonal blocks K^^^ (S) = 1, K^^^ (5), ■ ■ -, K^-^^ (S). Then the action of Sp (2, 3?) in A has 
the concise description 



lj^Kj{S)tj, t' = K^'\S)k^'\ (6) 

We are now in a position to present the generalized uncertainty principle. For each J = 0, 
|, 1, ■ ■ ■ form the square matrix Clj, of order ( J + 1) (2 J + 1), with operator entries, through 

the definition (^j is a row vector with the same entries as the column vector ^j) 



(7) 



We may write Clj in more detail in the block form 

/ 



1 r' 



It is to be understood that each element of j is written as a linear combination of the 
Tjn,nS using @). For purpose of illustration, we detail one of these blocks: 



3,0 



2,1 — -^-L 1,0 



T' I ih rp 

1,2 + "2"-'0,l 



Tr 



0,3 



Let Mj = (Vt) be the hermitian c-number matrix obtained from Vtj by taking (entrywise) 
quantum mechanical expectation value in the given state p: 



Mj 



(8) 
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It will prove useful to write Mj in the block form 



/ 1 MO'i ■ ■ ■ M"'^ \ 
M5'0 M5'5 ... Afi-^ 



J,. J 



(9) 



.■'\t 



where M^'^' = (^"'^^^ ' ) is a {2j + 1) x (2/ + 1) dimensional block, and M^''^' = (m^-^ 

Since C = 1, = 1 for all states. For purpose of illustration, we write out a few leading 
blocks of Mj explicitly: it is clear that the row vectors M^'a and M*^'^ have entries {q,p) 
and {q^^qp^p^) respectively; further 



qp + 



ih 



r 



q^p + ihq \ 



q^p - qp"^ + f p 



\ qp^ — ihp p^ 
and finally, the 3x3 hermitian block M^'^ has the form 

/ q^p + ihq^ q^p^ + 2ifiqp — \ ^ 

q^p — ihq^ q^p^ + ^ qp^ + ihp^ 



Here, g'^p" = {Tm.n) stands for the average of q^p^ with the Wigner distribution as the 
weight as in (|^). In other words, Mj is the matrix formed out of the moments of the Wigner 
distribution function, of order utmost 2 J. 

We now prove the important fact that the nonnegativity of the density operator p forces 
the hermitian matrix Mj to be a nonnegative matrix, for every J. For a given fixed value 
of J consider the operator 

V = 5Z $Z '^i* ' 

j=0 s=-j 

where arbitrary c-number expansion coefficients which can be arranged into a 

(J + 1) (2J + 1) dimensional column vector C. Now form the operator 



(10) 
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which is hermitian nonnegative by construction. Since p > 0, we necessarily have tr(p > 0, 
for every choice of the coefficients {cj^s}- But from (|^), ( [T0| ) we find 

That is, tr(/5(^) = C^MC for every C. This completes the proof that p > implies Mj > 
for every J. 

A little refiection should convince the reader that this is the generalized form of the 
uncertainty principle we have been after, and we state it as follows: 

Generalized Uncertainty Principle:- Let Mj be the hermitian c-number matrix formed out of 
the moments of the Wigner distribution of a state p in accordance with the prescription (^. 
Then 

Mj>0, j = o,l,l,--- (11) 

For a given state not all moments will exist in general. It is clear that in such a case where 
Mj is finite only for all J < Jmax, our generalized uncertainty principle should be modified 
to read Mj > 0, J = 0, |, ■ ■ -, J^ax- 

While the hermiticity and unit trace properties of p are refiected in the reality and nor- 
malization of the Wigner distribution, the generalized uncertainty principle presented in 
the concise matrix form (|lT]) exhibits the constraints on the moments of the Wigner 

distribution resulting from the nonnegativity of p. While the conventional uncertainty prin- 
ciple is such a constraint on the first and second moments, ours is a generalization to all 
orders. It should be appreciated that the canonical commutation relation enters Mj in (pH]) 
through (H). 

The following mathematical lemma is helpful in analyzing the content of this generalized 
uncertainty principle: A hermitian matrix Q of the block form 

Q = 

is positive definite if and only if A and B — C A^^ are positive definite. The proof simply 
consists in recognizing the congruence 

Qr^Q' = LQL\ 

where Q' is a block diagonal matrix with diagonal blocks A and B — CA^^C'^ . 

The usual uncertainty principle (|l]) is contained in ([Tl|) as a particular case: it is equiv- 
alent to the condition detMri > 0. Next consider the case J = 1. Use of the lemma with 

( = {M^^i AfO'i) renders M,j ~ M^, where 





M'j 



/ 1 
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Now Mj > implies Mj > which in turn imphes that its diagonal block Ma '2 
M2'^ > 0. Written in terms of the moments, the last condition reads 



/ q2 _ q2 qp — qpj^^A 



\qp — qp — ^ — 



> , (12) 



which is precisely the usual uncertainty principle (|I]). 

One more application of the lemma on the nontrivial part of Mj further strengthens the 
positivity requirement on the other diagonal block M^'^ — M^'° M°'^ to 



C= (m^'^ -Mi'°M°'5) . (13) 



This 3x3 matrix condition, together with the 2x2 matrix condition (|T2|), constitutes a 
complete statement of the generalised uncertainty principle involving moments of all order 
upto and including the fourth. 

It is clear that yet another application of the lemma, starting with Mj^a, will lead to a 



positivity statement on a 4 x 4 matrix which, together with ([T^) and ([T3|), will constitute a 
complete statement of our uncertainty principle on moments of all orders upto and including 
the sixth. Evidently, this reduction algorithm based on the above lemma can be continued 
to any desired value of J, and hence upto any desired (even) order of the moments, eventually 
rendering Mj block diagonal. 

We see from (^, (^ that Mj transforms in the following manner under G Sp (2, 3?): 

S: Mj-^Kj{S)MjKj{Sf . (14) 

The nonnegativity of Mj is manifestly preserved under this transformation. Thus, our gen- 
eralized uncertainty principle is invariant under linear canonical transformations. Further, 
the reduction algorithm suggested by the lemma is invariant under linear canonical transfor- 
mations, for it follows from (^ and ( p!^ that M^'^ transforms to K^{S)M^'^ [Sy under 
5GSp(2,3ft). 

An evidently useful way of reading (^^ is that the components of Mj, just as the Tm,n^-i 
transform as tensors under Sp (2, 3?). And the fact that our generalized uncertainty principle 
is invariant under Sp (2, 3?) means that it is implicitly stated in terms of the invariants of 
these tensors. These invariants, in the classical case, have been studied in great detail by 



Dragt and coworkers [|10 



While the nonnegativity of p implies the nonnegativity of Mj for all J, it is of interest 
to know if nonnegativity of Mj for all J implies nonnegativity of p. Phrasing it somewhat 
differently, we may ask: Given a real normalized phase space distribution whose moments 
satisfy the condition Mj > 0, for all J, does it follow that the phase space distribution is a 
bonafide Wigner distribution? 

From the very construction of Mj, it is clear that tr(p(9) > when O is of the form ( 
in (pUp. By linearity, this is true also when (9 is a (convex) linear combination of operators 
of this type(with nonnegative coefficients). Thus, ([TT|) will be sufficient to characterise 
the Wigner distribution if the set of all such convex combinations is dense in the space of 
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nonnegative operators. Intuitively, this may appear to be the case. However, the monomials 
Tm,n are generically noncompact, and hence a careful analysis of the issue of convergence 
should be made before one can make any claim in this direction. 

We have already referred to the intensity with which current experimental research deal- 
ing with measurement of the Wigner distribution is being pursued Since measurements 
are always accompanied by errors of various origins, it will be of interest to see to what ex- 
tent the Wigner distribution reconstructed in a real experiment respects the generalised 
uncertainty principle. Further, it may be of interest to examine the possibility of incorpo- 
rating these fundamental inequalities in the algorithm for tomographically reconstructing 
the Wigner distribution from measured data, in such a way as to improve the reconstruction 
itself. Finally, our analysis applies equally well to any other quasiprobability, provided we 
choose suitably ordered momomials and modify the product formula (^) accordingly. 
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